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SUMMARY 


When a rocket is launched, there is a short initial period of ac- 
celeration during which the rocket is unstable. As the flight velocity 
increases, the aerodynamic forces acting on tie fins and stabilizers 
become large enough to give stability. Various methods have been em- 
ployed to stabilize the rocket during this launching period. Guide rails, 
“gero length" launchers, booster rockets which produce high initial ac- 
celeration, and auto-pilot controlled nozzles are typical devices that 
have been used. 

This is an investigation of the requirements of a nozzle control 
which would stabilize the rocket during the launching period. The con- 
figuration investigated is unique in that the nozzle of the rocket is 
mounted as a compound pendulum, and the movement of the pendulum 
is utilized to furnish the signal for the nozzle control servo-mechan- 
ism, thereby eliminating the need for gyroscopic elements in the con- 
trol system. The pendulum motion of the nozzle caused by a change in 
flight attitude of the rocket is introduced into a computer which pro- 
duces an output signal proportional to the attitude of the rocket. This 
attitude signal is fed back to the nozzle control, which positions the 
nozzle. 

The results of the analysis showed that the rocket was unstable 
during the take-off period when the nozzle control acted on the rocket 
attitude signal alone. Stability over a narrow range of feedback gains 
was indicated for the system using a simple lead circuit as a nozzle 


control, or in other words, when the nozzle control acted on both the 
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attitude signal and the rate of change of attitude of the rocket. The 
damping characteristics of this system were poor. By changing the 
nozzle control function to include a response to the acceleration of the 
rocket attitude, the damping charactcristics were improved and the 


range of feedback gains was widened. 
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I. INTRODUCTION 


The rocket stability problem has existed for many centuries. 
Historians record that rockets were first used in battle by the Chincse 
in 1232 A.D. Sy 1400, rockets were in widespread use as major wea- 
pons by the warring nations of Europe. The poor stability of the rocket 
in flight was responsible for its gradual replacement by the cannon 
which gave more accurate trajectories. In the late 1700's, India em- 
ployed rockets with improved stability characteristics against the 
Dritish. When Dritish soldiers attempted to defeat the forces of Hydar 
Ali, Prince of Mysore, a trained corpse of rocket bombardiers launched 
rocket missiles at the charging Sritish cavalry, inflicting many casual- 
ties. The Indian rockets were made of iron cylinders about eight inches 
long and two inches in diameter, with a bamboo shaft about ten feet long 
attached to give stability. 

nm English military officer, Colonel William Congreve, became 
interested in the Indian advances in rocket development, and worked 
to further improve the rocket for British usec. Congreve devised a 
launching system consisting of a tube mounted at an angle on a tripod 
base, similar to the mortar. Long shafts were attached to the rockets 
to keep therm flying in a straight linc. The British employed rockets 
against Napoleon at Boulogne, and in other bombardments with great 
success. 

During the War of 1812, “.mericans were first introduced to 
rocket warfare when the Hritish bombarded Bladenburg, Maryland, and 


subsequently Ft. McHenry at Baltimore. An American by the name of 
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William Hale improved the stability of the Congreve rocket by replac- 
ing the wooden shaft with three curved fins on the rear part of the 
rocket body. The exhaust gases were deflected by the fins, causing 
the rocket to spin in flight. These rockets were used in the war with 
Mexico, and were found to have superior trajectories to the stick sta- 
bilized types. Again the rocket declined in popularity partly because 
of its inferior trajectory compared to that of the gun. No other out- 
standing developments in rocket stabilization were made following 
Hale's work until early in thie century. 

With the advent of the airplane, and an increase in knowledge in 
the field of aerodynamics, the stability problem of the rocket in Night 
became defined more clearly. The use of tail fine, or tail fins with a 
shroud ring as on the Bazooka, and later, the use of controlled airfoils 
solved the in-flight stability problem to a great extent. 

In the first part of the rocket's flight path there is a destabilizing 
effect from the change in mass of the rocket as fuel is burned. The 
acrodynamic forces must counter this tendency in order to stabilize 
the rocket. During the launching period, the aerodynamic forces are 
too small to counteract the effect of the change in mass, and the rocket 
ie unstable or at best neutrally stable. Consequently, stability during 
the critical take-off period must be attained by some means other than 
the use of aerodynamic surfaces. 

Various means have been employed to stabilize the rocket during 
the launching period. Guide rails to hold the rocket to a straight flight 
path until it has gained sufficient speed to become aerodynamically 


atable have been used successfully. The length of the guide rail has 
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been reduced to a minimum in some applications by firing the rocket 
with a very high initial acceleration. The German V-2 rocket was 
stabilized during the launch by an automatic pilot which controlled 
four carbon vanes in the jet stream behind the nozzle. The American 
Viking" and the German "Enzian'' employed a controllable nozzle 
positioned by an automatic pilot. 

A unique means of launching a large rocket with a booster rocket 
mounted free to swing as a compound pendulum was investigated by 
J. C. Norris in 1951 following a proposal by H. S. Tsien of Cal Tech. 
(Ref. 1). For the system investigated, the booster rocket was mounted 
below the main rocket on gimbals off the center of gravity of the 
booster as shown in Fig. 1. A deviation of the main rocket from its 
course would cause motion of the booster relative to the main rocket 
causing a change in the direction of the thrust axis. It was hoped that 
a practical configuration could be devised to give launching stability. 
The results of the investigation were somewhat discouraging in that 
for a typical vehicle-booster combination for which stability was indi- 
cated, the required length of the tail boom was too great to be prac- 
tical. 

Since the basic idea of the pendulum mounted thrust unit appcars 
sound, the possibility of utilizing the pendulum in conjunction with a 
servo-control is suggested. 

In this thesis, the requirements of such a servo-control system 
which would stabilize the rocket during the launching period are inves- 
tigated. The rocket has a controllable nozzle mounted as a compound 


pendulum; conceivably, the motor and nozzle of a liquid propellant 
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rocket could be mounted as a unit as shown in Fig. 2. The motion of 
the pendulum with respect to the rocket is fed to a computer which cal- 
culates the attitude of the rocket. The attitude signal from the com- 
puter is fed to a nozzle control. The aim of this investigation is to 
determine the requirements of the nozzle control such that it will sta- 
bilize the rocket during the launching period. 

The equations of motion of the configuration were first derived, 
assuming all aerodynamic forces to be negligibly small during the 
launching period. Motion in one plane only was considered. External 
forces considered were the rocket thrust, assumed to be constant; the 
jet damping force arising from rotation of the jet stream; and the gra- 
vity force. The effect of the change in mass of the rocket resulting 
from the burning of fuel was considered. The displacements of the 
rocket and nozzle about the initial condition were assumed to be small 
quantities for the period of the stability investigation. 

Following the derivation, the equations involving the pitching 
motion of the rocket and the motion of the nozzle were transformed to 
simple algebraic equations by use of the Laplace Transform. From 
the transformed equations, it was possible to solve for the transfer 
functions of all of the components of the control system, as shown in 
Fig. 4, except that of the nozzle control. The form of the nozzle con- 
trol transfer function was estimated, and the stability of the resulting 
system was investigated by use of the Nyquist Stability Criterion. The 
Evans Root Locus plot was used to find the damping coefficients of the 
modes of oscillation of the system. Successive estimates of the form 


of the nozzle control function were made in order to improve the 
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stability of the system. 

It was found that for a nozzle control function employing feedback 
of the rocket attitude signal only, the system was unstable. For a con- 
trol function employing feedback of both the rocket attitude and rate of 
change of attitude, it was determined that a stable range existed be- 
tween upper and lower limits of feedback gain; the limits depending on 
the constants of the control function. The damping of the system was 
poor. 

The nozzle control function was revised by introducing additional 
lead circuits, or sensitivity to the accclerative change of attitude of 
the rocket. The damping was improved by this change, and the range 


of gains was widened. 
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Fig. 3. Block Diagram of a Feedback Control System 


A typical feedback control system which might be employed to 
guide a rocket of the controllable nozzle type is shown in Fig. 3 above. 
The attitude of the rocket, as indicated by the angle 9s which the 
rocket makes with the vertical, is measured by an attitude gyro. The 
gyro signal is matched with the desired rocket attitude signal and the 
error, or difference, is fed to a nozzle control. The nozzle control 
applies a moment // to the nozzle which results in a relative angular 
displacement Oy 9s of the nozzle with respect to the rocket. The 
thrust axis rotates with the nozzle, causing a change in the rocket 
attitude. 

Now consider the case of a rocket with the motor and nozzle 
mounted as a compound pendulum as shown in Figs. 1 and 2. With any 
change in rocket attitude, the pendulum moves in a calculable manner. 


From a knowledge of the system dynamics and atime history of the 
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pendulum motion, the attitude of the rocket may be computed. Suppose, 
then, that the angular displacement of the nozzle from the center line 
of the rocket, Gy -%s » is measured and fed as an electrical or me- 
chanical signal into a computer with an output representing the attitude 
of the rocket gz . The gyroscopic element is thereby eliminated 


from the circuit. 


F2z F Os out 
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Fig. 4. Block Diagram of a Feedback Control System 


Employing a Compound Pendulum Nozzle 


The block diagram of the proposed feedback control circuit as 
shown in Fig. 4 above may be interpreted as follows. 

The desired input of the coordinate 2s : , or angle which the 
rocket makes with the vertical, is matched with a feedback signal rep- 
resenting the actual angle at any given time. The difference between 
the desired angle (in this case it is desired that Gs =O for vertical 
takeoff) and the computed angle is fed into the nozzle control mechan- 
ism as an error signal, ig - Is . The error signal initiates the 


application of a moment /77 ‘to the nozzle gimbals. The moment 
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contributes to the nozzle dynamics to yield a certain position of the 
nozzle relative to the rocket, which is represented by 94 — 43 . The 
thrust line is rotated off the center line of the rocket by movement of 
the nozzle; this contributes a moment and side force to the overall 
rocket dynamics, resulting in a final rocket position. 

The functions Fe . Fs , and Fa appearing in the block diagram 
of rig. 4 are transfer functions derived from the equations of motion 
of the configuration which have been transformed by the use of the 
Laplace Transformation. The transfer functions are linear operators 
and can be defined by inspection of 7 block diagram of the system 

e., fe = 4s ; {= Fh O,-0; ). The nozzle control trans- 
fer function fF, is the unknown - a which is the subject of this in- 
veatipation. 

The e uation relating the output of the system to the input may 
be derived as follows: 


Following the circuit from Q, iy (60 an OUT 


Q, ouT [Qs in Q; |F, re 


.round the loop 


Ny 


ae “10,4- Qy FFF, 


Q - F Fe Fi Qs in 
. | + FR Fe 


Substituting for Q. in tho first equation 


Qsz out QurpuT — I - Fi Fa Fa Fi Fs 


Qa in INPUT ia FR (OER TPR 
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The System Transfer Function therefore becomes: 


K e = OuTpuT _ Fi fh Fs 
ss ~ INPUT | + Fife Fe 


The stability analysis thus involves 1) derivation of the equations 
of motion for the rocket and movable nozzle configuration, 2) conver- 
sion of these equations by the Laplace Transform to solve for the trans- 
fer functions fF, , Fs , and Fi , and finally, 3) synthesis of a nozzle 
control function F , and investigation of the complete system trans- 
fer function |. G, for its stability characteristics. 

The Nyyuist Stability Criterion and the Evans Root Locus Method 
are applied in the stability investigation. Refs. (2) and (3) contain ex- 
planations of the use of these and other methods of analysis applicable 


in the field of stability and control. 
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Ul. EQUATIONS OF MOTION 


The assumed rocket configuration with all external forces to be 
considered is shown in Fig. 5. A schematic diagram of the coordinate 
system and dimensions of the configuration are shown in Fig. 6. 

Frictional and aerodynamic forces are neglected. The change of 
maes of the main rocket body with burning of fuel is considered. Mo- 
tion in one plane only is considered. 

It is convenient to employ coordinates 4, through “ie in the 
derivation, but of these, only four coordinates, q, through Va! are 
independent. The coordinates 4, and q. measure the horizontal and 
vertical distances respectively to the center of gravity of the entire 
configuration, and have the dimensions of length. Coordinates +; and 

4, <imstant the angular displacement of the main rocket center line 
and the nozzle line respectively from the vertical, and have the dimen- 
sions of the radian. 

The cquations of motion in the four degrees of freedom, namely 
in the directions of coordinates ;., » y 


2 
using the generalized Lagrangian momentum equation: 


4 (32 }- st . 4. (1) 
atl éq, _ 


Elimination of Extraneous Coordinates 


b Ge’ and t6 are derived 


By taking moments about the center of gravity of the system as 
shown echematically in Fig. 7, the relation between the system center 
of gravity and the centers of gravity of the main rocket body and the 


nozzle may be found: 
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Fig. 7. Sketch Showing the Relation Between the Centers 


of Mass, M, and M2, and the Center of Gravity of the System 


FMe.g= Mg (4.- 4,) — Meg (4,-4,)= 0 


M | (9° q.) = M2 (¢.2 47) (2) 


From the geometry of the system as shown in Fig. 8 the follow- 


ing relations can be obtained: 
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Fig. 8. Geometry of the Configuration 


Equating the vertical distances between mass centers: 
%e 7 ke 4, 8 (‘h. §) 


_ 4-4 
Ne = Ces (q- §) (3) 


I 


E\juating the horizontal distances between mass centers: 


™ V5 = 4, Sin (9, - §) 


S 2 ts 7 47 
+ Sin(q4,-S) 


By considering similar triangles in Fig. 8 it follows that: 


(4) 


oe 2. St te 


ee 


es , +, - 44 
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Mi(¢e-%2) _ Ma (42 4e) 


5 

mt: ie - 4.) M. (4,-4.7) 
Substituting Ey. (2) into (5), the result is 

Mm. (Ge? ty = Ma (4,- 48) (6) 


Eqs. (2), (3), (4), and (6) are four linear expressions containing 
the four extraneous coordinates ts j Ve? t7 ; tLe and the principle 
coordinates 4%, 4,., +3 and Va . The former may be expressed 
in terms of the principle coordinates as follows: 


From Eq. (3) 











1.= 440s (9,-§) + Ye (7) 
From Eq. (4) 
t7= ts - Ag Sin “a0 S) (8) 
Substituting (8) into (2) 
ts= 4, t+ a Sin (4° - §) (9) 
Substituting (7) into (6) 
Mm, 
Ve= Fa japan, 14 6°8(45" §) oo 
Substituting (9) into (8) 
+7 + ¥ ai 7 §) (11) 
Substituting (10) into (7) 
16 = tat -§) (12) 





It will be noted that the above equations contain the variables 


Ag and § from Fig. 7 which must be expressed in terms of the 





47 = 


principal coordinates. 


In Fig. 8 by the law of cosines 


QR. = VRE + Ay - 2 Bg coo (yp, - Oy) (13) 


and by the law of sines 


y _ Q 4 
i> i Sin (44-4, ) 


Sin, < sin(3,-45)| 


1} 


8 


-| XM, 
- _ = ' 5 —- Sth = 14 
hy 6e= & Oe lz stm (4 i) “? 
To simplify Eqs. (13) and (14) consider the following approximations: 
Pid Sin +; Cos ¥. = | 
Then from Eq. (13) 


ha = ~J/f,* 14-28 a =~ -% 


From Eq. (14) 


fr 


ts ra +; 7 (4, - +3) 
ve Q Wi Hi 
13795 sa "¥ +4 
and since 
ke = &, -A42 
| Xo | 
4- - $254 a an 7 i. = — - in(q,- 8) 


The approximations then are 
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Xe 
Sin (q,- s)=A eyle Q, -Q. +4 


cos (4,-8) = (15) 
Ri — he, -~Ay 





Substituting (15) into (9), (12), (11), and (10) yields the following re- 


lations between the extrancous and the principal coordinates: 











¢.- 1, Wie (2.43 - 42 4,) al 
4e= 4, 4ek (A- 4) 7) 
Y= brie, (4 95- f4e) = 
1 —_ ~42) a 


The Kinetic Energy Term 
The kinetic energy | appearing in the Lagrangian equation of 


‘ motion may be derived as the sum of kinetic energies of translation 


and rotation by referring to Fig. 6: 


=41,4, +4 1,4, +¢™. (454%) thm (45 +4) 20 


Derivation of the Equation of Motion for the Coordinate ., 
The first term of the Lagrangian equation may be obtained by 


first taking the derivative of Eq. (20) with respect to +, 


oT » 345 e A 
—_- — M , + M ne 
Sea in! TY 
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The extraneous coordinates Ys and %, appearing in the above expres- 
gion may be eliminated by differentiating Eqs. (16) and (18) and substi- 


tuting the results. Thus we obtain: 


Pa (mM, +Mz) a - MM (£:4,-42 





385, Mth 
The above expression may be differentiated with respect to time to 
yield: 
4 (27) (m4m,)§, + mg, — Mam (gg ge \y MMe lg g | 
dt er )= owe. 1 Mi+M, ( 4; 4) am (4 iz Q, in| Mah 


The second term of the Lagrangian eyuation may be obtained by differ- 
entiating Ey. (20) with respect to 
2” = 0 (22) 
t, 
The generalized force on , can be computed by considering the 
work done by the external forces as the coordinate + is varied, all 
other coordinates remaining fixed. 


In gencral, 
Work; = Aki dG, 5 of kt = =a! 
L 
In particular, for coordinate +, , the force ae becomes, referring 
to Fig. 6: 
g - Work _ (Fin dy + Deora) O41 
.= 


“¢. 43) 


With the approximations of Eq. (15), the above expression becomes 


2 ,= Fq, + D (23) 
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Substituting Eqs. (21), (22), and (23) into Ey. (1), the equation 


of motion inthe 4, direction becomes, finally: 





(Mit Ma) 4, +M, 4, — Seer (4i43-44,) + pean y (4.45- 94) 


~ Fa, + D (24) 


aé 


1 $ 
Derivation of the Equation of Motion for the Coordinate ¢. 


In a similar manner, differentiating Eq. (20): 


ot : 04. ; 24e 
—— = Mm, > ‘ + M £ 
4, a) a 24, 


and, differentiating Eq. (17): 


re = * ae 2 (4, - - £3) 


(M.+M M,+M2)" 

04. 

04, 
Similarly, differentiating Eq. (19) 


+= 41- Ma M (4 f,) 


(mi+™Ma)- 


= | 





846 = I 
3 +, 
Combining the above results: 


= (mM +™2)4, - Mam (4-42) 


mM, + M2 





rs 4. 
By differentiating this expression with respect to time the first term 


of the Lagrangian equation becomes: 
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d fal) _ , oe fg Ma Me " 
rt ite = (M.+Mz2) Set, &, rr: | (4, §, ) (25) 


The second term of the Lagrangian equation is: 


a = (26) 
te 
The generalized force ot, , referring to Fig. 5, may be ex- 


pressed as: 


2 - Work, — [-(M.+M.)g t F cosy, - Dsing, | 44, 
-_ ————— — ee —————=— 
a4. 44, 


With the approximations of Ey. (15) 
2 = -tmom)g + F-O4, (27) 


Substituting Eqs. (25), (26), and (27) into (1), the equation of motion 


in the +. direction becomes: 





(m.+™M2) 3, + M, 4, pe (fh) = F -9(™,+™,) + bq, (28) 


\' 
Derivation of the E juation of Motion for the Coordinate 43 


The first term of the Lagrangian equation is obtained in an iden- 


tical manner to that employed for coordinates 4 and Ye : 


Differentiating Eq. (20): 








oT - 0a 
ei. _- + mds se +5 + Me V7 —t! (29) 


: +3 la 
and, differentiating ~ (16): 

15 +, + 7 M,tMa (Xf, + fea.) - 

04. _ Mi Q, 


04, M, tM 


oe A 


ip (Ys- 4,4) 
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also, by differentiating og. (13): 


I= 4.7 the (ha,- Re) et (hg ey) 





M. tM, 
04, = ee M. e 
O45 = oo te 


Substituting into (29) 
ay. 
ot ~~ 1 a M, M2 0, q - MM, 0, fh, q 
oF, , PA, My, 3 Mit M, # 
Differentiating the above expression with respect to time, the first 


term of the Lagrangian cquation becomes: 


™ i ee MoM. * mM, (, 
eS) =! 43 nat Wot (4, ds- 4 ie) 


Ma MM, 
* Garam? 43 Fe) (30)" 
From Eq. (20) 
i, (5,28 i 65 ) + M, (i, 3% rt. ~~ 2 te me (31) 
43 94 3 
By differentiating Eq. (16): 


Vs = q, * oe th. 43- M, +4) - permits eH) 


te = ™ Mi mM,£, (32) 
2¢ « - (M.+Ma)° 





Hy differentiating Eq. (18): 


 F- ” fees ts- ®4,) Eesha M44) 








29, _ MM, (33) 
24, ~  heeeaenaey” 


ie a, ) oe 

The term 147 ’ which accounts for the propellent mass tranéfer out 
of the rocket is not correct, inasmuch as this moment should be calcu- 
lated from the velocity with which the propellent leaves the transfer 
tubes and enters the motor. Corresponding crrors appear in subse- 
quent equations. The numerical effect of the resulting terms is negli- 
gible. 
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, a " 
In order to obtain aM, and = , it becomes neces- 
3 





3 
sary to differentiate the exact expressions for 5 and te: since the 
4; term has vanished in the linearization of the exact expressions 


qs. (10) and (12) will be differentiated and the results then linearized 
From Eq. (12) 











M2 a 
to= tothe, *4 == | (© d 4,-§) (12) 
qe . 7 ~ ae Xe sine 3 - fast a Ay cose —_— Q, cos © 
84, - . Bu. 21s,,.434- §, AaB 6 - ae 2. 
oq ~ MtM, 34, tatM, dG, M,tMa 3%, 
3 
. Ma M, 2X4 30 e. Mam, {, d cose Mi Bs 6 r=) 
(M, tM)" 9 OV 4 (Mit Mi.) tM.) r) + 3 mM,tM, a +3 
* Ma g 0 cose (34) 
M,+Mez 043 
From Eq. (13): 
2, = L* +h, ~2R, £, 605(44 -4,) 
By differentiation: 
OX4 X, Re Sin (44-43) 
— + ea. (35) 
+3 Q 4 
ake = f = A, Aa [sin(4,-45)] | 44 -4s] (36} 
Dat 4 —$ 
4 
+ 6 «a Sek (9 -4.) cos x 
DE Oh) Gil wed tee) 
f, 
ae -4 Ka 








(37) 
1 (44-43) 
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From Eq. (14): 
dae | 
= ", — Sin [4 Cia (44-43) | 
By differentiation: 
rs co$ (Q,- Ga). See- 43) — (2nd order Lerm) 





b= 4, - 
3 Ji- Oto he) 
5] - Ga, 
So %; Le (Fy 43) _ 
6 tds 
Sap 7 Ge (Ho WG Fs) = 
go heel a= (Seen) Red} 
34, i- [BSG] Le 
068 ww X2 ~ XZ, 
ay ~ + 22 ay 
843 7 Re a - 
Also: 
OSinoe _ oo ~ ay 41 
‘Tl prot Ep x4 - 
o cos 6 -s Sines ~ ANG hg) (42) 
04. Fy By 


By substituting Eqs. (35) through (42) into (34) and ignoring terms of 


second order or higher in the smal) quantities: 














Ove Mi 0 ¢ f, ' 
~- 43) 
MtM, ae M, — = $3 
34, a 
In a similar manner, from Eq. (10): 
p= & = Sine o ¢ cose O44 _ ML R4cos MMe ¢ cos 
Vo m Ve — . MitMy ot MAM. . ‘ vy ‘ .)* . 


By differentiating this expression: 
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ay Mi 24 sine 28 + Qe 6 dsine + Sine 6 Sl 
24, ~ MM, 34, 043 




















OY 3 
_ My ts [coro - ~ Re sine 22 = M, My 5 O- sine lat ® 
MitMa a °Fy) Mam. | oq, 845 
MM, aX4 PY. 
4 —*4 cose — Sino 44 
(M+ M,)* Ae "3 35, | (44) 
Again substituting Eqs. (35) through (42): 
0%5 M, M, M, ? 
= —_—_—__— 45 
oq, MitMa 4 +3 * (ma M,)? ‘3 lat 
Having evaluated o%e and Ste from the exact expressions, 
$3 ys 
at can now be obtained by substituting these values in Eq. (31). 
Y 3 
7 no __ M, ‘ mM. mM, 
24, MyM . ** (mM, +m) 3 (4.9,-% t) 
_ Mem 9 (9-9 4 MM fi 
Gamage 4, eg, 7, (46) 


The generalized force 


a a Work 3_ 
7 * &3y 


Then, referring to Fig. 6 


= ci a 
= “ign 4 mig ste + F sin gy 5 tr (Lark) sing, 
3 3 


+ Fees, 3 [Berli ria)eosga] + Deng, (Get) sng] 
x 3 
= D sit 1439, ie (hr h)cos 4, 
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Substituting in values of os and te from Eqs. (li) and (10) 
2 : 2 | 
2: [fein Desa] [a.-fadeeee the) 


i | Fees Mei D Sin +4 | om ie mea pi (4 +4) cos9,] 





Performing the differentiation indicated: 


= ‘ _ mM Sine aft 
) nae “oe of 
t | Feos4, — D sin 4+ | i. ae 395 “<7 J mad 


Now, substituting Eqs. (35), (41), and (42): 


= = (Fq, +0) {- fe 1h . Say hos a4) (- 14. ows))) 
+ (F 44) | hoy f : ba £4 ale $s uA Gas] 


Finally, retaining only firet order terms, the expression reduces to: 


43> - we LF (44-45) +] wily 


The equation of motion for the coordinate 15 can then be ob- 








tained by substituting Eqs. (30), (46), and (47) into Eq. (1). 


oe M, M, 2 2 =! 
\, be ° Rh (hi 1 f.4,) + rs Pk a (44° fh +4) 











MiMi = um4, yet ; 
+] 4s "Mar be a. ts he imam)! (ha, 44) 
mim £1 (2-42) _ ives i 
. (M.+ Ma) * by = 7 Mam IF (44 ,)+ ‘ou 


| 





eels 


‘a “a 


Derivation of the Equation of Motion for the Coordinate 44 
By differentiating Eq. (20) with resocct to: +, 

er 
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From Eq. (16): 


Ys = 4,+ Gtk, (Ligrtg,) “fe (Fi 93° M44) 


a . . 34 
= Jig + ™ Is 3a, + Mi +75 (49) 








~ifferentiating: 
; = v 
o 44 7-7 


From Ey. (18): 


1, = 4.7 ae dy 4. 4,)- feof bs & 4,) 








Differentiating: 
3 ‘ _ M, Q, 
a Qa cai 
844 “— e 


Lubstituting these juantities into Eq. (49) 


ol = = mM, M, 4 ” » Pe ’ 
o44 = 4 * MitM, (84, +4) 


The first term of the Lagrange equation is obtained by differentiation 


of the above expression: 
mM, mM, 22 ee oe MoM, fe = . 
(35 =1, tte r MitMae (4.45 42 4,4) 5, (M.tM,)” (44; 144) sate 
From Eq. (20) 
My hr pg O48 
BF, = Gs + 4G) + “(He tia) 60 
Py differentiating Eq. (16): 


ry * me (4 19° 4,4) me (4 ts 4s) 
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0g cs , MLM, Qa 
OV — (mM +M, )* 
And by differentiating Eq. (18): 
; M2 _ Fe M, 
49 = 2, + gm, (419g 24 ) (M,+™:)” (4.9,- 414,) 


04, 2 Mimi 4a 
ode x (Mt Ma + 





S46 and ore 
044 oy 4 


in the case of the Ys coordinate, to differentiate the exact expressions, 


In order to evaluate it is necessary, as 
and then make simplifying approximations: 
From Eq. (12): 


MoM, 


+= :..- rT) ek Sin O a {4 Cose +t f, Cos 6 




















Differentiating: 
9%. Mr ok SinO 6 - =. hq g2ine g — 2 NG sive SE. 
2%, MM OY 4 Ts ++ Mirth " 
- =». oX¢ C638 pa gcone mh Se Oke cos 8 
(MM) 4 (M.+Ma) o%4 Mth 844 
Mi 0 0c0S@ 
+ 
Mam Ba, (52) 


By differentiating Eq. (13): 
ke _ fda (44-43) 
9%, Le 


Also by differentiating Eq. (36): 


bie  2F,2: (G4 - 4s) 


= 


44 RL 4 








By differentiating Eq. (14): 

. . . = 
ts X4 
And by differentiating Ey. (38): 


38 = QQ, (44-43)(44 43) 
we 3 


Ce 


O/ 


94, X4 
Then 9% may be evaluated by the substitution of the above ex- 


844 
pressions and Eqs. (35) through (42) into Eq. (52): 





OF _ As 4.4, - MiM.£2 


344. ‘tM (M.ae)* +4 - 





From Eq. (10): 


a M, oe a lh aXe 
te = qi Yes. Se sin6G nam 6% (M.tm)" 





Differentiating: 
Pe) — Mi Dy sino Z® + 046 922, sine odes _ st ee 
ie +4 0% 4 + tC 8%, 





rere _ Mi M, dcos 6 re) o 2k ' 
ob, | (atm) f ne. > (54) 





Then cas may be evaluated in a simila manner to 2% by 
a o+ 4 
substituting for the extraneous terms in Eq. (54): 


o%¢ a mM, 42 oe Mr M42 
54, i mtMe V4 (mM. +a) * 44 — 


The second term of the Lagrangian equation may be evaluated 
by substituting the derivatives obtained in the foregoing steps into 


Eq. (51): 
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oT ag M, mM, b Ma m, £2 , 
oq, wah YO Ram Vt 
_ MiMi d mi mite 
(mitmi)$ wee (hs h 4, * OM 3 (2. 4) 4, (56) 
The generalized force 
- Wor K a 
Ay = “ay, 


By reference to Fig. 6, it is evident that: 


As = (F si Yat Deest4) 35 {9 ,- (fa +P) sim | 


+ (Fcosq, — Dang) 2, [Spe ~ (fa rFs) cos 4, ] 


3 re) 
wigs — mia gs 





Substituting for the extraneous coordinates and simplifying, as in the 
derivation of Eq. (47): 


1M 
2,2 - > (ph 6) " 


Introducing Eqs. (50), (51), amd (57) into Eq. (1), the equation of mo- 


tion for coordinate 4 becomes, finally: 
vs at Mi Mi fa mM £2 
Tn 447 aan (hhh) eae (4a 8 je) 


a Mi Mi £2 : Mim, 4a me 2 
M. tM a -1y © mM tMa be 4. ta (M.+M,)> ay (£195 - h 1+) 


: mnt (8-4) 4, = - on +4 s] i 
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IV. COMPUTATION OF SYSTEM COMPONENT 


TRANSFER FUNCTIONS 


Only the two equations involving the rotation of the rocket (co- 
ordinate +; ) and the motion of the nozzle (coordinate V4 ) are con- 
sidered since the stability investigation is made for rotation only. 
First, the thrust, jet damping, and moment of inertia terms are sub- 
stituted in terms of the equation variables and constants. 


Rocket thrust: 


F= —Mq_1sp 
Jet damping force: 
LD - eds +4 


(The negative sign is necessary in order to make the term 
positive, since M, is a negative quantity. ) 


Moment of Inertia derivative: 


a. 
i = Fin R, (R = radius of gyration) 
.= = wR + 2h, 
, a ; 
Iz MR, (since A, will be small for 


most configurations) 


Substituting the above terms in Eq. (48): 
Mm, mM, & 


i, o. as MiG Z, + ™+Ma ene’ (F3 4, )+ mets (Ag ‘Fs h ds) 


- ems, (l.q5-%q4) + pa Q-f)g, + Mam’ 5. 


MitMa 


™ . ' ' 4, 
-ea Fee > =e [- Migtsp (34° 4)7 Mids +] (59) 





% z 








a52. 





Again, substituting for a ; D , and qT , and introducing the applied 


* 
moment into Eq. (58): 


om 1 A. iM, La : ; 
Ig EEK (gehen) — A (5,44) 





_MaiM 4 (4,5 ,- fh. 44) _ Mim, Mi. Li (6. - ~h)g,- MiMi 9, g 


MitM, (M,tM.r)> TM), 
tM f Yahi = mM. 4s 4, (fy + Me & )+ (60) 


The Laplace Transformation is applied to the equations prior to 
the stability analysis. Ref. (2) explains in detail the use of this trans- 
formation in such analyses. In essence, the transformation converts 
the differential equations to simple algebraic equations. 

Since the vertical and horizontal velocities q. and q, are 
quite emall] during the launch, the terms involving +, and 1.4, are 
neglected in the equations. 


Application of the Laplace Transformation to Eq. (59) yields: 


MitMe (Mit M) MitM. (MtM)3 


—_ mM, Ma 4s & * Mim Ae mim 44s 
= Oa) FARR )S + eee ys 


Qs} [1+ ee] s+ [inks My Mobs > 5 + [taint ew 


, fmumtigtse _ mem, 4h a 
mM, t Ma (M: +M,) (61) 


where Q; represents the transformed coordinate q; , and S§ is the 


complex variable. Similarly, transforming Eq. (60): 


os 

The effect of the control moment M on the main rocket is neglected 
in Eq. (59). For a system in which the nozzle is of comparable size 
and weight to the main rocket, the term should be included. 





233. 
os {uti s+ PAAR] s ~ ERREAKT) = Moot 


: rm hh 
Q, if. + =e — | a , [pasa mM fs (f, xt l \s 4 aa (62) 


where M represents the transformed moment ‘7. 








In order to solve for the transfer functions F, and F, it ie 


neceseary to rearrange the equations slightly; Eq. (61) becomes 


Q, \f,* mM, £, (8 ams + [t, >, 1 a, 4 (2, oma MM HAs Mm, 4. ks | 5} 


M, to. (M,+m,)* Ph .tMe 


M,My 2,2 mame he 4 Mi My hfs MM Agdse Mim hte 
= |Q,- Qs) fea a |S + Gacy * M,t Ma ar MitM. = (Mem | ata 
Similarly, Eq. (62) becomes: 


Q, lerccscra ) a)s4 pout a. Is (h riers) s} 


MitMr (mM, +m.) * Mit Ma 


= oy~ Oy) ft BE] + [ght intel) 


M, Mm, Sf, g, 
1 peas a aaa a 


Then the transfer function Fs from Eq. (63) becomes 











[rat S (BS mb fe + mM, &) s+ [tei er Ms i. “+ 


(Mite) Mit Ma 


Qs 





ete a | 
Q4- mM 0,( 2-20) - [ ie mitt, £(0-4)_ mM. 4 Ls 
Sn OO lela, =| ver | 


By substituting iq. (63) into (64), the function Fe may be written as 
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= Q4- Qs : | 

OS ea a ey, 
a ee — —— — Fe eS eS SS SS a nee 
a fz + Mena 


Pails at 0. - MM Lg (66) 
fal, tina |S me + zm f, (fs +8 + m)} S (Mi, Q)3 
As an aid in the numerical computations for the stability analy- 


sis, the transfer functions are symbolized as follows: 


F, = i — -£ 


= —_ = (67) 
Q4- Qa sg 
> O1-Gs .. | 


;- xX \ _ (68) 
Mm (Xp -W 





where: 


was 
| 


Gi: esc 4, b * mimo £5 S + mm, 4, I, a M, moh 
G =| + jeu (Mit Ma)? Mt Mi Mit Ma (mM, th)3 


“. 





z, 


een Meth. (™M,+Ma)* ate 


monte h) x 2] 5° [mi ma(t 6) ow) (e+e) sf 


m4, 0-4) rm fi(f- 4 AAs 
= {fie ne re maatal |S + [ae + mam ES nds hts} 


ie SS mh) +(e vfs (Be BR) 7? saat] 
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V. SYNTHESIS OF THE NOZZLE CONTROL TRANSFER 


FUNCTION "Fy,", .ND STADILITY INVESTIGATICN 


To illustrate the determination of the nozzle control transfer 
function, a specific numerical example will be worked out using the 
German V-2 as e2 typical vehicle. 

“8 indicated in the Analysis, the system transfer function may 


be written as 


4,Gs5 = _ & Fg (69) 
1® & Fy %& 


“.lso0, the computer synthesizes the dynamics of the rocket to give 
fy = fy (70) 
Substituting Eq. (70) into (69) 


a FA 
<4"? 4.7. (71) 
lt #, Fa Fs 


Let 
(72) 


Substituting for -, and fF, from Eqs. (67) and (68): 


AG - ae 73 
& bette / (73) 


se a numerical example, consider the V-2 rocket with pivoted 
rocket motor and nozzle. The parameters y , » , 7 , andW 


have the following values as calculated in the .:.ppendix: 








=86-. 


ne 2675 G&G - 46/5 = 4 Of' nce 


7= G2 emo §”- //45 


2@ £86030 8% - /465 85 (74) 


me €8#2 6* € 10.858 


The problem now is to determine a function fF, which will give 


the system suitable stability charactcristice during the launch. The 


procedure will be to: 


1. 


é.. 


i.esume a function fF, 

Check the system for stability by the Nyquist criterion. 
Having established a stable system, determine the roots of 
the system transfer function by means of the Evans Root 
Locus Method. 

Knowing the roots, at a typical operating gain of the feedback 
system, investigate the damping of the system. 

If the resulting motion is unsatisfactory, improve the choice 


of fF, and repeat the process. 


In explanation of the methods to be used to investigate the sta- 


bility of the feedback control system, it will be noted that the equations 


of motion are linear differential equations with constant cocfficients. 


The solution for the particular coordinate in question can then be rep- 


resented in the general form 


§,t Set 
%; ia a e +C,e eS eas 
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where s isa simple root of the equation. In order to have a damped 
motion which will subside with time, it is obvious that the complex 
roots ‘S; ’ must all have negative real parts. -..nother way of stating 
this requirement of a stable system is that all the roots ‘Ss; ” of the sys- 
tem must appear in the left or negative half of the complex plane. 
The Nyquist Criterion 

The Nyquist Criterion is an application to the simple feedback 
network of the Cauchy Theorem in the theory of complex variables. 
This theorem may be employed to show that if “Cay {s the vector rep- 
resenting the function of a complex variable “s, the function having 
‘™” poles and” zeros within a closed contour C in the ‘s” plane, that 
as the point S moves around the contour once in a clockwise direction, 
the vector | CO» carries out 7-7 clockwise revolutions about the ori- 
gin. If the contour C te chosen to include the entire right half plano, 
and the function in question is plotted as the complex variable encircles 
the right half plane, the difference between the number of zeros and 
poles of the function that exist in the right half plane will be indicated 
by the number of encirclements of the origin. Either the number of 
zeros, or poles existing within the contour must be known previously. 


In applying the theorem to the feedback network, the function in 


question is usually expressed as 


fe AG Cs} 
—_——s« - -—_ °°  — | 
[+ kG) 


where the function _ Ge. can be factored; thereby determining the 


zeros of the function y 4 G." . The denominator “/+ #6" in general 


cannot be factored so readily, but application of the Cauchy Theorem 





=-~§8< 


will show the number of roots of the denominator existing in the con- 
tour chosen, or the right half plane. 

In order to obtain a plot independent of the gain ot, the inverse 
of the function may be considered. 


kG, *G 


‘ 
Then if the vector 4s Gs encircles the origin, the vector = 


will encircle the point “-/' as the variable “s’ moves around the con- 
tour of the right half plane. Or similarly the vector e will en- 
circle the point as jl The function Zs } is then plotted as ‘Ss en- 
circles the right half plane. By application of this Nyquist Criterion, 
the stability or instability of the system can readily be determined. 

If the nozzle control transfer function s " is first assumed to be 
equal to a constant (F, 5 K,) the system transfer function takes the 
form 

fe Gg= tS. 4h Va) = 
(+ 4 ya yw) 

For the Nyquist diagram, the inverse of GC is plotted as the vec- 

tor S encircles the right half of the S plane. The function to be 


plotted then is 


ne F822 Si 
Ors) x 


In order to use the Nyquist stability criterion it is necessary to know 
whether the function 4G has any zeros in the right half plane. By 


factoring Eq. (73) the following expression {s obtained: 
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AC) = ti [- 0 OY), (Ei Cae (77) 
s/s -. 006 f{s- (~. 0204+ 6.25/[5-(. Oz 04-¢6.28)/ 


Since the function 4~G'‘has one zero in the positive real half of the 
plane, the Nyquist Criterion requires that the plot of die must 
encircle the point - 4 once in a counter-clockwise direction as the 

§ vector circles the right half of the s plane once in a clockwise 
direction. 

Fig. 9 is the Nyquist diagram of Eq. (76). As the frequency 
increases from zero to positive infinity, thence to negative infinity and 
back to zero, there is one clockwise encirclement of all gains from 
O tooo. The net encirclement may be considered to be two clock- 
wise rotations, indicating two roots of the system transfer function 
existing in the right half plane. The system is unstabie for all gains. 

Next, assume the nozzle control transfer function to be of the 


form f= Ki (7+ 7S) - then 


ki Cr Tedd yet) 


i & = 
4k Ch BN ae) 


(78) 
The Nyquist diagram for the above cquation is shown in Figs. 10 and 
11. Again the number of zeros of “G “existing in the right half plane 
determines the Nyquist criterion. The factored equation now becomes 
; $- (a s4/979/(S-/97 6/ 
VEG 
Ky = Tt finer 0 MS a ta. 
(e.25)) 


/s- 006 |/s-f- 0204+ ozo ie 


Since there is one zero in the right half plane, again it is desired to 
4d ¢ 
have one counter-clockwise encirclement of the gain 4 as an indica- 


tion of stability. 
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It was found necessary to plot the Nyquist diagram to several 
scales to indicate clearly the crossing points of the negative real axis. 
The enlarged portions were plotted for the case when 7=/_ _, and ap- 
pear as Figs. 12 and 13. From these plots it is clear that there is 
one counter-clockwise encirclement of all gains from “%</< to 

A= 75°. The configuration is stable for gains in this range. 
The Evans Root Locus Method 

This method is a graphical solution for the variation of the roots 
of the feecback system transfer function K 6 as the gain ‘K ie 
varied from zero to infinity, where 


ao Gn 
ent 2 stm 
l+ 466) 24s 





The plot is constructed by observing the following rules (Ref. 2): 
1. For 4:20 , the poles of ge are the poles of Gy) 
2. Yor fm aq _, the poles of a ral are the zeros of C6) ; 
with any missing poles taken at s-q@. 

3. The product of the distances from a root to the poles of AG : 
divided by the product of the distances from the root to the 
seros of 46 ts equal to the gain “¢ :' 

4. The eum of the anglcs from the zeros of AC “to a root, less 


the sum of the angles from the polcs of Li 3 “to the root must 


equal 7 , or "77 ''plus an integer multiplied by 97. 


The locus of the roots of hy Gs " can be determined graphically 
by trial and error, observing the above rules. The gain scale can be 


plotted along the path of the roots, and for a chosen operating gain, 
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the roots can be determined by inspection. 

The root locus plot of Eq. (18) for the case when 7=/ is shown 
in Fig. 14. The roote emerge from the poles of *G at zero gain and 
proceed along the loci indicated to the zeros of *G ata gain of infin- 
ity. Gains at intermediate points can be approximated by dividing the 
product of the distances to the poles of 7G by the product of the dis- 
tances to the zeros, and by the constant gain factor in the equation 


(49.2x 107°), It is apparent from the root locus plot that the range 


for a stable configuration is from the gain where the roots “S, % 
cross into the stable half plane to the gain where the roots he 
cross into the unstable half plane, referring to Fig. 14. The loci in 
this range are plotted to an enlarged scale and appear in Figs. 15 and 
16. As an example of the computation, approximate expressions for 
the gain along the loci in this range are: 


For root 


7 E.25- w)(E.25%4,) ley 2) 


2 
0 22050 ‘call 
,) CY9 2x10 D0 9.79N09. 64) V7 Fas se eet) 


For root 


Laz 02 0. 25NC2. sKé.25:) (for a ern 


4:,)? ————————— 
WEG VEF. [0 525 ig /+ FH) 92.2 y/o 


hej FETOW+ O2 ov) 
It will be noted that the gains at the limits of the stable range 
check within limite of accuracy of the scale used on the graph as shown 


in the table below for the Ny uist and the Evans Root Locus plots. 
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Upper Limit of K Lower Limit of K 


Nyquist plot 75.0 
74.80 


Having established the gain scale along the root loci in the range 










Evans Root Locus 


of positive stability, referring to Fige. 15 and 16, it is possible to 

choose an operating gain and determine the values of the roots Sx 
te Vl 

and ry, at that gain. Some possible operating points have been 


tabulated below 


-.0031 +i. -.00325 +i. -.0033 +i. 
-.0031 - 14.156 -.00325 -i .171 -.0033 - i .185 


-.0068 + 16.25 -.004 +16.25 ~.0003 +1 6.25 


-.0068 - 1 6.25 -.004 - i 6.25 -.0003 - i 6.25 





The low frequency roots S, = are associated with the motion 
of the main rocket, whereas the high frequency roots = a are as- 
sociated with the motion of the nozzle. The time to damp to half am- 
plitude and the period may be computed for the modes quite simply; 
the results of such a computation for a gain of 4=6O are tabulated 


below. 
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iZime to damp to 
half amplitude 210 seconds 171 seconds 


~eriod of one 
cycle 


36.38 seconds k second 





ee 





Number of cycles 


to half amplitude 5.7 cycles 171 cycles 





It is apparent that the damping is quite poor, and an improve- 
ment is desirable. The acceleration of the V-2 at take-off is roughly 
equal to 30 ft. jpew.”. If it is assumed that at a velocity of 200 ft. /sec. 
the aerodynamic forces stabilize the rocket, the period of take-off in- 
stability is asproximately 7 seconds long. Yor the foregoing example, 
the rocket conceivably could complete only one-fifth of a cycle of 
slightly damped motion in the take-off pcriod ae a result of a unit dis- 
turbance at the instant of firing. The nozzle oscillations are so rapid 
as to conceivably cause overloading of the control device. 

In order to improve the damping characteristics of the system, 
referring to the root locus plot, Fig. 14, the addition of one or more 
zeros along the negative real axis would serve to pull the roots 5, 
and 5. into the negative half plane more quickly, since zeros attract 
the roots much as lines of force are drawn to magnetic poles. This 
should increase the magnitude of the negative real part of the roots at 
a given gain. Also, the angle of departure of the root S, would be 


increascd counter-clockwise as a result of the angle condition of the 
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root locus plot. The root would then be forced to travel a longer path 
before passing into the unstable half plane. 
The addition of zeros in the left half plane involves adding addi- 
tional lead networks, which in real components have transfer functions 
C fa 7s) ‘ P 
of the type + €TS) where € + > 7G =, 
Assuming a new nozzle control transfer function to be of the form 
EAT SHOOT SN +75) 
—— r z 
(1+ 355 S) big s) 


4¢ 4% 
where /7=/ , tho factored expression for 4G is now 


46 = —o/ peal MEP 1-9 7008.64 _—_ 
S(S=006//5 PLOKS 7 ofs-( OL OV re 6. 2uf[steo2 Oy-< bz) 


The root locus plot of Eq. (33) is shown in Fig. 17. The range 


(83) 


of gain for etability is now from K = 4 to K = 918. 
By approximating the loci of the roots in the stable range by 
straight lines, the roots associated with various gains can be estimated 


conservatively. Some representative values are tabulated below. 











-.02 +1 15.25 
-.02 - 115.25 


If the damping time and period are computed for the latter 
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system, the following results are obtained. 


K = 600 Rocket Oscillation Nozzle Oscillation 


Time to damp to 








half amplitude 15.3 seconds 34.2 seconds 
Period of one 
cycle 11.6 seconds -412 seconds 
Number of cycles 1. hy dens 83 cycles 


to half amplitude 





Interpreting the results as before, the rocket now completes ap- 
proximately two-thirds of a cycle of oscillation during the launching 
period, but the effect of a unit disturbance is damped to a much greater 
extent. The nozzle oscillations occur at a higher frequency but the 


damping has been improved considerably. 
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VI. CONCLUDING REMARKS 


The analysis shows the feasability of stabilizing a rocket during 
the launching period by utilizing the motion of a compound pendulum 
suspended nozzle as the input to a feedback servo-control system. It 
was determined that the nozzle control for such a system must include 
response to both the attitude and the rate of change of attitude of the 
rocket, as computed from the pendulum motion, in order to obtain 
stability. It was further determined that changing the nozzle control 
to include response to accelerative motions of the rocket improved the 


damping characteristics of the system. 
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APPENDIX 


The V-2 Rocket 


















( 

: | 3500 Ky , saso Ks 

2. Ky | \ 1 

‘ { 

1 

EXPLOSIVE , 2s : + @ : + + + 
AVXKILIARIES ALCOHOL : OXYGEN | PUMPS 
t 


CG 
30.2 
37.5 
ai? 
OVERALL HEIGHT SO 
Weight and Balance Table 
Moment 
Wt. in Kg. Wt.in Lbs. Moment Arm, ft. lb. ft. 
1. Explosive 930 2160 7.42 16, 000 
2. Fuselage 1750 3860 26.65 103, 000 
3. Pump 450 993 37.5 37, 200 
4. Motor 550 1214 41.7 50, 600 
5. Auxiliaries 300 662 14.0 9,300 
6. Alcohol 3500 7720 20.83 161, 000 
7. Oxygen 5250 11600 30.2 250, 000 
Totals 12, 980 28, 620 727, 100 
Moment arm of C.G. = 124219 . a5 4 @. 
28,620 


Note: Data was obtained from Ref. (4) and is only approximately 
correct. 
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Calculation of the Moment of Inertia 





j a oe 
ashe 45 i 
Pi 

we 

A 

| ow 

[ om 

| W) 
| 
| 
! 
ub 2 & 


Assume the rocket is a solid cylinder as shown in the above sketch. 


Radius of gyration of a solid cylinder about its C.G. is 
[ae ke 
Rs 2 12 
(TWwac)* + €83.35)* 2 
BeOS Jaa = VIL a 


Additional Data from Ref. (4) 


F =Total Thrust = 27.2 tons = 54, 400 lbs. 


W = Total Weight = 12,980 Kg. = 28, 620 lbs. 


M, = Mass Flow =— 19, 320 Ibs. _ 193.2 lb. /sec. = —6 Ilb.sec. /ft. 
100 sec. 
: _ F_ . 54,400 . 
Ip = Specific Impulse S Mig = T93.2 = 282 sec. 
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Calculations for the Constants of Section Ul 


2 
ME 2? 406 ibs - os) SS 
2 
Boe ieit be: 37.7 oS 
i £ de . _-¢ lbosec 
M, = 193.2 3B = —6 MBAS 
k = 97 ft 


k= 2 Ft (radius of gyration of motor) 
, = 20 

}, 5.7 $F 

Rh = 0.3 ft 

I, = 4,0 ft 

I =M™Mk = 80,600 Ib ft sec” 
TI, = Mk = IS! Ib ff sec 
Tgg = 262 Fees 


: 
M,+M, > 806.705" 





li : ° 2 
4 = MM, bbe e* +- IA, & % MM, GA S MM, Aglsp mM, 44) 
+ roe, m, + M,) M+M, 





28 620 Ibs 


Xx = Jensee de 7 thsesy (20 HG 7##) 











" 2 
ibs) (_¢ 5 ME) (0G £1)(3. 741) _ Gs tae I tb oe "a —6.0 Pats M2 oft) (oft) 


(868 ? lb sec* )2 886 7 > sect 
Fr fr 


aan > Je pec & 
er 


y,* 2 
(37.7 b22* ) (-6.0 FE) (2oftM 3.7 44) 
(885.7 ees")? 
+r 


x = | 2675 S — GGI S— 1044 xX 10% | (16 ft) 


r< 


Mit, 8, (U-fer) 62 - 32 MP (RRA) 1 ae 
z 2 Lay t zu = Mt, M, . 
G bi Mr My ro (" K, (A,+ My)? v3 )s 


(888. 7 = - 


Y = (co.e0s lbft se , (65 Bee )(37 7 Wey og iL s* } 


(62a “\(o.7¢4) + 7.7 BF) (-6.0 Naot Meo —3.7)6t $706 


ue a lb a i 





(as, B38 ye epee ts 


888.7 lbsec* 
ft 





ef = | 92,380 5° — No | 
Zz 


MM, i (4,-) _ 2 MoM, ’,(£,-£,) 
\ M+ hr, I.) ~~ ( (M, +Mz)* +H As As tet A}}s 
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“_ Pwerne 20-37) _ 1) ay (22 (6320-37) | 


886.7 (888 7)* 


(37.7)(3 7) 

(-6)(4)( 34 3867, js | 
7 = | (2161-151) S* + (—652--11) 8 ] [6 ft 
vA = [ 2030 s = 8 S ] 

mM 2. 2 mM. M, t , M 
WW ~y &, + wate ) id (A = fs (4 + a, “))s 7" 
Ze 
MM, #4, 
(M,+M,)3 | 
- (as') (37.7) (3. 7) (37) (-6)(3.2)" 422. 7)(3.7) 

W = ( *— a,7 ~~ )s* +( (86a.7)t €)(4)(.3 886.7 Is) s 


_ (32 7)" (-6) (20) (3.7) 
(888.7)! 


VW = (si +495) s* + (-1I5 +11) S —.0054 


\A/ aw, [e+e s + 10.85 | 
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